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ABSTRACT 


Variation of the alcroiMve Intensity and spectrin due to gyro-synchrotron 
radiation from senl-relativistic particles InJ^ted at the top of a closed 
magnetic loop has been described. Using the recent high spatial resolution x-ray 
observations from the HXIS experiment of the SMM and from observations by the 
VIA, It Is shown tluit the high microwave brightness observed at the top of the 
flare loop can come about If 1) the magnetic field from top to footpoints of 
the loop does not Increase very rapidly, and 11) the accelerated particles 
Injected In the loop have a nearly Isotropic pitch angle distribution. The 
limits on the rate of Increase of the magnetic field and/or the average pitch 
angle depend on the geometry and location of the loop on the solar disk. 



!. INTRODUCTION 


Close correlation between the observed tenporal variation of the hard x*rays 
and microwave radiation during the impulsive phase of solar flares indicates 
that the same (or a closely related) population of electrons are responsible for 
both of these radiations. It is assumed that a wide energy spectrum of electrons 
is produced during the in^ulsive phase* with the lower energy ones (e < 100 keV) 
being responsible for the bulk of the observed hard x-rays and the higher energy 
ones (c > few hundred keV) producing the microwave radiation. However* one of 
puzzling results obtained from the SMH HXIS experiment and the high resolution 

ground-base microwave observations has been that the x-rays (15 to 30 keV) 
and the microwave radiation (at wavelengths of 2 to 6 cm) are not coming from 
the same region. In the few events observed, the tendency is for the bulk of 

x-rays to come from the foot points of flare loops (Hoyng et al 1981) 
while the microwaves are observed to come primarily from the top of the 

loop (Harsh and Hurford 1980; Harsh et al 1980; Kundu et al 1981). 

The purpose of this letter is to investigate the limitation that these 
observations impose on the . parameters of proposed models of 
flares. 

In general, x-rays which are produced by bremsstrahlung are simpler to 
analyze and give more direct information about the characteristics of the 
accelerated electrons. In almost all models the x-ray intensity is expected to 
increase from top to the foot points of the loop because of higher densities 
at lower regions (see Emsiie 1981). 

On the other hand, the microwave radiation, produced by gyro-synchrotron 
process, depends on the pitch angle distribution of the electrons, on the magnetic 
field structure and is affected by various absorption processes (Ramaty and 
Petrosian 1972). In two recent works we have developed a simple description 
of the variation along a magnetic loop of the pitch angle and energy distribu- 
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tion of electrons Injected at the top of the loop (Leach and Petrosian 1^1 » 
hereafter referred to as LPI) and have derived simple analytic formilas for 
evaluation of the gyro-synchrotron flux from $em1-relat1v1st1c particles with 
arbitrary pitch angle and energy distributions (Petrosian 1981, PI). We use 
these results to evaluate the variation of the microwave flux along a single 
flaring loop and compare It with observations. In this letter we consider the 
general features of the models. A more detailed analysis will be published 
elsewhere. 


2^ nnam(MW«iW^ 
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II. GENERAL DESCRIPTION 


To simplify our analysis we consider a flaring loop of width D» length wR 
and a central magnetic field line which Is a semicircle (see Figure la). We 
assume that the strength of the magnetic field B Increases from top to the 
footpoints (symmetrically), B « B^b(s), so that the width also varies as 
[? “ D^/b. This means that other (non-central) field lines will deviate from 
the assumed circular configuration but If D « ds/dInB. such deviations can 
be Ignored. If the distribution fj(u»Y) In pitch angle cosine y and 
energy y nf electrons along the loop Is known (fd^dyds Is the nimd>er density 
of electrons within y to y + dy, y to y + dy and s to s + ds), then 
the gyro-synchrotron emissivlty j^(s,e) at frequency v and In the direction 
e with respect to the field lines can be calculated using the equations In PI; 

j^(s.e) = (e^VjjSlne/c)(v/VjjSine)** fj(3^cos0,y^) G{0.y^) . (1) 

Here, G(0,y ) Is a slowly varying function of 0 and, in general, like f , 

decreases with Increasing y^, Vj^ = eB/2mnc = 2.8 x 10® Hz(B/gauss) Is the 

gyro-frequency and y^ (the energy of the electrons with the highest contribution 

-2 

to the emissivlty, B* ” 1 - y^ ) depends on the energy spectrum of the electrons. 

For electrons with pitch angle distributions which are not extremely anisotropic 

-6 

and have a power law energy spectrum, f 5 (yiy) « (e^.- 1+ y) .or for a thermal 
electron gas with tmnperature kT (In units of mc^), f(y,y) « y(y*-l)** 
and to a good approximation (cf. PI); 


|4v/[3V|j(« + l)s1n0], s 1 

l(2vkT/vjj)(l + 4.5vkTs1nVvjj)“’^^, kT < 1 


( 2 ) 
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In the optically thin regime, v > v*, the observed Intensity Is 

I^(s,e) « Jj^(s.0)dl » j^(s.8)0^/b\lne. v > v* (3) 

As frequency decreases, rises till the critical frequency v*, where the 
optical depth is of order unity; then with further decrease In v (optically 
thick regime), It begins to decrease. As mentioned above. It Is not clear 
which absorption mechanism will be the dominant one. Except for self-absorption 
the other mechanisms depend on the ambient plasma paran^ters. Because of model 
dependence of these other absorption processes, we will consider only the self- 
absorption process, In which case the Intensity can be approximated by 

1^-2 mv*, for v < v* . (4) 

Here ^E^ is some average particle kinetic energy (in units of mc^, E = y - 1). 

For a Maxwellian distribution ^e} = kT and « v*. For a power law distribution, at 
v»Vb*(^E) « Yq - 1 s V** so that The critical frequency v* is a 

complicated function of v^, 9 and D; v* = v^sinS H(Vb»B,D). However, in most 
cases H is a slowly varying function of the parameters. We will Ignore its 
variation and set v* = H v.sine. 

0 D 

As is evident from eqs. (1) to (4), the spectrum and flux of the gyro- 
synchrotron radiation depends primarily on the particle distribution function f. 

Thus, before we can calculate the emissivity we need to specify the variation 
of density, pitch angle distribution and energy spectrum of the electrons along 
the field lines. This Is the most complex part of the problem and the part 
which depends strongly on the assumptions of the models, on the many unknowns 
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of the flare p1a»na and the acceleration of the particles. Below we shall 
consider a few forms for this distribution covering most of the proposed models 
for the Impulsive phase of solar flares. 
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III. SOME MODELS 

The characteristics of the models depend on the distribution in phase 
space of the accelerated particles and on the variation of the magnetic field 
and plasma density along the loop. 

These properties depend on the combination of average pitch angle a^, the 

quantity RdlnB/ds and on the dimensionless column depth dT^ “ l4Trr^lnA(E+l )/E^Jnds 

“13 

(n is the ambient electron density, r^ = 2.8 x 10 cm and InA « 20). 

For example, if the electrons, injected at the top of the loop, have small 
pitch angles (a^ « 1) and the magnetic field is nearly uniform (RdlnB/ds « 1) 
a beamed thick target model will be the result. The x-rays of energy E then 
will be produced primarily at the regions where the column depth measured 
from the top of the loop, is of order unity. This is because the number of 
electrons with kinetic energy E decreases rapidly when exceeds unity 

(cf. figures 3 and 7 of LPI). The HXIS observation that the 20 kcV x-rays 
originate primarily from the footpoints indicates that t < 1 throughout 
most of the loop, which means electrons are injected at a column depth 
N = ^nds < lO^^cm"^ above the footpoint (or the transition region). Note 
that for high energy electrons needed for the microwave radiation this means 

Te < 0.001 . 

In the other extreme case, if particle distribution is nearly isotropic 
(a^ of order unity) and the magnetic field varies rapidly (RdlnB/ds » 1 for 
« 1)» then the electrons will be trapped and radiate x-rays primarily from 
regions with the highest magnetic field and ambient density which again will 
be the lower parts of the loop. In general, the parameter which determines the 
degree of beaming or trapping is 



sin^a^ ^InB/dt^ . 


(5) 
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For Cg « 1 one has a beamed model and for » 1 a trapped model . 

In the low density regions of the loop (t^ « 1)» ?£ » 1 and the dlstrlbu- 
tion of particles are determined by the adiabatic Invariance of B/sIn a. For 
example, for an Injected spectrum at x « 0 of ^^(sln o,E), the distribution 
along the loop becc*L*s (cf. LPI, eq. 7) 

fj(sin^a.E) = fQ(s1n^ot/b.E) (6) 


On the other hand. If the density is large so that x^ exceeds unity 
much before the magnetic field has changed significantly, then the particle 
distribution is determined by the collisions with the ambient plasma. In 
general, there is ro simple analytic expression for the distribution except 
In the small pitch angle regime. For example, for Injected electrons with 
gausslan pitch angle distribution and energy spectrum f^lE), eq. (18) of LPI 
for relativistic energies gives 


%(p.E) 


f„(Ec) 


2aexp I- 




Inc 


(a* + Inc), C “ 1 + T, 


(7) 


For non-relativistic energies Fq(Fc) ■* Fq(Ec )/C • 

We now consider the microwave radiation frcxn a few models with different 
val ues of the parameter c^* 

1) Uniform Trap Model . Let us first consider the simplest (but scxnewhat 
unrealistic) model with nearly uniform magnetic field and isotropic pitch 
angle distribution. Such a model will result if x^ « 1 and if particles 
are injected isotropically at the top, so that according to equation (6) 
fj.{y,L) * course, to have a trap model, the magnetic field 
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must eventually vary. We assume B = const throughout except at the foot- 
points where B -»■ « rapidly. 

Since che magnetic field is uniform ® const), then, according to our 
earlier discussion, v* and (y* - 1)*^ will vary as sin9 so that the spectrum 
and intensity of the microwave radiation will depend primarily on sin0 and the 
frequency of observation. In Figure 2a we show schematically the variation of 
V*, and with cos9 at a given frequency. At higher (lower) frequen- 

cies, the Yq curve is shifted to higher (lower) values. The quantity ^j(Yq)G then 
changes accordingly, decreasing the increasing y^. In Figures 2b and Zc 
we show the expected microwave spectra at a few values of 9 (9 * tt/ 2 for the 
highest curve). For power law spectra, the spectra cross each other because 
in the optically thick region “ Y^ “ i » while for a thermal source of 
uniform ttuperature, I^ « kT. 

As evident from these figures the variation of intensity with e will 
depend on the frequency. As shown in Figure 2d, at v > the maximum 

turnover frequency, the optical depth « 1 for all 9 so that 
I a f X G (cf. eqs 1 and 3) decreases rapidly. For v < the intensity 

V % f r tMX 

increases slowly (as Y^"^ » power law) or remains constant (thermal) for 

jcose] < ^°^®crit '^*^®crit^ ' ^ ^ ^°^®crit 

< 1 and the intensity decreases rapidly again as f x G. 

The angle 0 and consequently the variation of the intensity along 
the loop depends on the angle n = s/R (cf. Figure la) and on the location 
and orientation of the loop on the sun. In general, for a loop near the 
solar equator, at heliocentric longitude i> 


±cos6 = cos^sinn + cosies in(t>cosn , 


( 8 ) 
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where \i> Is the angle between vhe projection of the loop and the solar equator 
(see Fig. lb). In general, the location along the loop (i.e., the value of n) where 
the intensity if maximum (cose » 0) will vary with ♦ and ip. 

For a loop near the center of the solar disk (^ % 0 and all \p) cos8 * ±simi 
so that the intensity will be highest at the top of the loop. As we move awa> 
from the center 4> > 0 , the maximum intensity (8 = ir/2) occurs at f 0 

(e.g. n * ♦ for r|» * 0). However, when projected on the solar disk the 
maximum intensity will appear approximately half way betweeen the footpoint'J. 

^ote that for ♦ f* 0 the microwave emission will noc necessarily be symmetric. 

This kind of configuration rather than an asymnietric field geometry, may be the 
explanation of some observed asymnetries (Alissandrakis and Kundu 1978; Kundu 
and Vlahos 1979). 

This picture will change near the limb (4> = tt/ 2) where for ^ = tt/2, 6 « tt/ 2 and 
the microwave intensity is uniform all along the loop or for =0, cos8 ® cosn 
and the highest intensity occurs at the footpoints. However, the observations under 
consideration here with two distinct footpoints do not refer to these configurationi., 

We conclude, therefore, that ihis model agrees with the microwave observations. 
The basic reason is that it is the value of the component of the magnetic field 
perpendicular to the line of sight which determines the brightness of a synchrotron 
source. Furthermore, as shown in Figure 2d, the variation of the intensity (along 
the loop) is strongest at the highest frequency, so that the source of the micro- 
wave radiation will appear smaller (more concentrated toward the middle of the 
footpoints) at higher frequencies. Comparison of the 2 cm (Marsh and Hurford 1980) 
observations and at 6 an (Kundu et al 1931; agrees with this aspect of the model. 

The flux from the whole loop can be obtained by integration of the intensity 
along the loop. As shown by the dashed line in Figure 2b or 2t at v > the 

total flux will be decreasing as in eq. (3) with most of the contribution 
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coming from 8 » it/2. At lower frequencies the contribution from other perts 
begin to become significant, and the total spectrum will be somewhat flatter 
than ( but will approach asympthotically to) the spectrum in eq. (4). This 
type of spectrum is commonly observed (Hurford, Marsh and Zirin 1981; see also 
Solar Geophysi al Data) and is attributed to inhomogeneities in the source or to 
other absorption processes (Ramaty and Petrosian 1972). We can see here that 
even a uniform loop (uniform in field st-ength and particle distribution) can 
qualitatively reproduce such observed microwave spectra. 

Note that the above picture remains qualitatively the same even for non- 
uniform magnetic field as long as the magnetic field variation is slower than 
1/sinO. Any such variation, however, would give rise to a slower variation of 
the intensity along the loop and to stronger asymmetries for loops away fm 
the center of the disk. 

2) Non-uniform Trap Model (Cg >> 1). Now we consider a nx>del where the 
magnetic field varies rapidly throughout the loop (RdlnB/ds » 1) instead of 
the extreme variation concentrated ut the footpoints of the above model. As 
mentioned above, particles responsible for tne microwave cmissioi. have << 1 
so that for of order unity » 1. As in model (1) the pitch angle 
distribution is governed by eq. (6), according to which the distribution 
broadens as the field strength increases. 

This model is different from piodel (1) not only in its allowance for 
non-isotropic distribution but, more importantly, because of its rapid variation 
of the magnetic field (and v^^), both v* and increase from top to lower 

parts of the loop reversing the trends of model (1). The dashed line on 

Figure 2a shows the schematic variation of v* (the decrease near cosO - >1 is 
due to the unrealistic circular shape of the assumed loop). This increase in v* 
(and the decrease in «nd the increase in f x G) is stronger for loops 
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away from the center of the disk. In this model then the top curves on 
Figures 2b and 2c will correspond to the footpoints of the loop such that the 
trend in Figure 2d is reversed, as shown by the dashed line which clearl/ 
disagrees with observations. 

3) Uniform Beamed Model (c^ « 1). Now we complicate the models by 
injecting electrons non-isotropically with small pitch angles. If the magnetic 
field is constant or varies slowly, then « 1- The variation of v* and 
in this model will be similar to that of model (1). However, the quantity 
f(3^cos0,Y^)6 will now vary not only because of variation of Yq but also 
because of the non-isotropic nature of the distribution and the variation of 
3^cos0 along the loop. For this model the variation of the distribution is 

_3 

given by eq. (7). As mentioned above, along most of the loop < 10 

so that if is greater than a few degrees, thenineq. (7) Inc < , 

a = cos”^(B^tor0) 0 and f^(B^cos0,E^) « exp(-0^/a^). 

For all orientations and locations of the loop we have at the footpoints 
7t/2 < 0 < tt and at the top 0 < 0 < tt/2 so that the value of f is, in 
general, much larger at the top, which makes the variation of 1^(9) in 
Figu.e 2d even steeper than that of model 1. There are, however, two problems 
with this result. First of all, the stronger the beaming the fewer the number 
of particles contributing to the flux. If the self-absorption was not important, 
this would be a natural explanation of the old discrepancy between the number 
of electrons needed for x-ray and microwave emission. However, with strong self- 
otion this model will produce less microwave flux than observed. The 
secoMi difficulty with this model is that eq. (7) is valid for small pitch 
angles while we are interested in values of f, typically, at ot :: tt/2. As 
inspection of Figure 1 of LPI will show, the number of particles at such large 
pitch angles can vary by large factors with slight changes in the value of a^, a 
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or in the field geometry. Thus, it is difficult to make a definite statement 
about the validity of this model. 

4) Non-uniform deamed Model (C£>>1). Finally, we consider a model 

diametrically opposite of model 1 in that the field varies strongly and particle 

Ditch angle distribution is non-isotropic. Since the magnetic field varies 

strongly, we have a situation similar to model 2 where v* increases with cos0 

(see dashed line in Figure 2a) and decreases with cos9. Here, f(3^cos6,Y^) 

increases with cos0 not only because of the decrease of (as in model 2) 

but also because of variation of 3 cosO. As evident from eq. (6), the pitch 

0 

distribution broadens with increasing field strength, giving rise to a larger 

value of f(P^cosO,Y^) at the footpoints as compared to the top. For example, 

for a gaussian injected spectrum in a loop at thi center of the disk 

_2 

f (P cosO,Y ) « exp(-a )« 1 at the top of the loop, while at the footpoints 
0=0 o o o 

where P cos6 = -1 and b » 1, f ,o is about unity. Thus, in contradic- 

O 0 =Tl/t 

tion with observations, varies even faster than that given by the 

dashed line in Figure 2d. 
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IV. SUMMARY 

We have presented a qualitative description of microwave emission from 
electrons injected at the top of a closed loop with particular attention to 
the variation of intensity and spectrum along the loop. As evident from the 
discussion of previous sections, many parameters enter into this description 
of the models. We emphasize here the effects of the orientation, location and 
geometry of the loops and the pitch angle distribution of the accelerated 
electrons. 

We have considered the total intensity (disregarding the polarization) 
for four models which qualitatively agree with the HXIS observation of 
'v^20 keV x-rays. We find that the high resolution microwave observations can 
be reproduced by models where the magnetic field increases slowly from the top of 
the loop to its footpoints at the transition region. Faster field variations 
give stronger emission at the footpoints. In addition, another requirement is 

that the accelerated particle distribution should be nearly isotropic. Model 1 
satisfies both these requirements. If the accelerated particles are strongly 
beamed along the field lines, then the footpoints will be brighter than the top for 
a rapidly increasing magnetic field, but for a uniform field the situation is 
uncertain and a more detailed analysis of the pitch angle distribution is 
needed. 

We have considered models with the extreme value of the critical parameter 

For intermediate values of this parameter results intermediate to those 
described will be obtained. 

We have neglected absorption process other than the self-absorption. These 
other processes will be more important at the lower, high density, regions or 
of the loop end, therefore, could reduce the intensity of the footpoints 
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(see e.g. Ramaty 1969, for effects of Razin-Tsytovich suppression 
mechanism). These and consequences of ou.' other simplifying assumptions, in 
particular, setting v* = H^v^^sInS , along with the expected variation of the 
polarization, will be described elsewhere. 
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FIGURE CAPTIONS 

Figure 1. Geometry of loops on the sun. 

a) Assumed geometry for a loop. Electrons are injected at the top 

n = = 0, 0. The indicated values of cosG is for a loop 

at the center of the disk. 

b) Loops at various solar longitude (J> and orientation For 

loops at the center of disk cose = tsinn, at the limb 0 = it/ 2 
throughout if = ir/2 and cos6 = cosn if = 0. 

Note that for loops with distinct footpoints the point where 
line of sight is perpendicular to the field line will appear midway 
between the footpoints. 

Figure 2. Schematic representations of variation of various quantities with 
angle 0 and frequency v for model 1. 

a) Variation of the critical frequency v*, energy and 6f(y^), the 
quantity determining the emissivity in the optically thin region, with 
cose (dashed line for model 2). 

b) Synchrotron spectra at different values of 0 for a power law 

elect.'^on spectrum. |cos0| = 0 (0 = t/ 2) for the curve with 

V* -V * and increase gradually for the lower curves. Dashed line 
max 

for the spectrum integrated over all angles 0. 

c) Same as (b) except for a Maxwellian electron distribution. 

d) Variation of intensity with cos0 at two different frequencies (dashed 
line for model 2). The location of the top and footpoints are shown 
for a loop at the center of the solar disk. For loops with = n/2 the 
footpoints move in the direction of the arrow as the loop moves from the 
center to the limb. At the limb all points along the loop are at cosO = 0. 



Vahe Petrosian 


Institute for Plasma Research 
Stanford University 
Stanford, CA 94305 
(415) 497-1435 



cos 9 - 0 


I 




(b) LOOPS ON THE SUN 


Figure 1 





